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It is shown that the massless energy spectrum of electrons and holes in graphene leads to the 
strongly non-linear electromagnetic response of this system. We predict that the graphene layer, 
irradiated by electromagnetic waves, emits radiation at higher frequency harmonics and can work 
as a frequency multiplier. The operating frequency of the graphene frequency multiplier can lie in 
a broad range from microwaves to the infrared. 



In the past two years a great deal of attention has been attracted by a recently discovered, new two-dimensional (2D) 
electronic system - graphene, built out of a single monolayer of carbon atoms with a honeycomb 2D crystal structure 
[H, [|J . The band structure of the charge carriers in this system consists of six Dirac cones at the corners of the hexagon- 
shaped Brillouin zone [3j,|4(, with the massless, linear electron/hole dispersion. The massless electron spectrum leads 
to unusual transport and electrod yna mic properties, which have been intensi vely studied in the literature, see e.g. 
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Electrodynamic properties of graphene have been theoretically studied in Refs^iiliilillil, [lOJSJJ^JiljMjii, 

lEEEEIliEEiEil!. The frequency dependent conductivity[Il,[iJ[2EJp,[23], as well as plasmon [2^ tH^, 
plasmon-polariton 24] . and transverse electromagnetic wave spectra 3l| have been investigated. In all these papers 
electrodynamic response of the system has been studied within the linear response theory (for instance, using the 
Kubo formalism, or the random phase approximation, or the self-consistent- field approach). In this Letter we show 
that, apart from all the fascinating and non-trivial properties of graphene predicted and observed so far, this material 
should also demonstrate strongly non-linear electrodynamic behavior. In particular, irradiation of the graphene sheet 
by a harmonic electromagnetic wave with the frequency fl should lead to the emission of the higher harmonics with 
the frequencies mfl, m = 3, 5, . . ., from the system. The operating frequency of such a frequency multiplier can vary 
from microwaves up to infrared, and the required ac electric field is rather low, especially at low carrier densities 
and low temperatures. The predicted non-linear electrodynamic properties of graphene may open up new exciting 
opportunities for building electronic and optoelectronic devices based on this material. 

To qualitatively demonstrate the non-linear behavior of graphene electrons consider a classical 2D particle with 

the charge — e and the energy spectrum e p = Vp = V^Jp x + Py in the external electric field E x (t) — Eo cos fit. Here 

V is the velocity of 2D electrons in the energy band (in graphene V ~ 10 8 cm/s 0, 0)- According to the classical 
equations of motion dp x /dt = —eE x (t) the momentum p x will then be equal to p x (t) = — (eEo/fl) sinf2t, and the 
velocity v x = de p /dp x is then v x (t) = —Vsgn(sm fit). If there are n s particles per unit area, the corresponding ac 
electric current 



jx(t) = en s Vsgn(sin fit) ~ en s V— { sin fit 
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contains all odd Fourier harmonics, with the amplitudes j m , m — 1,3,5..., falling down very slowly with the harmonics 
number, j m ~ \Jm. Notice that at the density n s = 6 • 10 12 cm -2 and at V ~ 10 s cm/s (parameters of Refs. [l|, |2j) 
the current amplitude j = en s V in our simple estimate gives a giant value of j — 100 A/cm. 

The above consideration does not take into account the Fermi distribution of charge carriers over the quantum states 
in the conduction and valence bands of graphene. To get a more accurate description of the non-linear phenomena 
in the considered system we use the kinetic Boltzmann theory, which allows one to get an exact response of the 
system not imposing any restrictions on the amplitude of the external ac electric field E(t). Using this quasi-classical 
approach we take into account the intra-band contribution to the ac electric current. The inter-band contribution to 
the electric current, due to the transitions between the hole and the electron bands, is ignored. This imposes certain 
restrictions on the frequency of radiation fl, which will be discussed below. 

Consider a 2D electron/hole gas with the energy spectrum e p ± = ±V^Jp x + p^ under the action of the field 

E = (E x ,0), where the sign + (— ) corresponds to the electron (hole) band, E x (t) = E§e at cos(flt), and a — > +0 
describes an adiabatic switching on of the electric field. Assume that the Fermi energy cf lies in the electron (or 
the hole) band and that the temperature is small as compared to ep, T <C ep- The momentum distribution function 
of electrons f p +(t) = f p (t) (we omit the sign + for brevity) in the collisionless approximation is described by the 
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Boltzmann equation 
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which has the exact solution 



fp(t) = F a {Px -po{t),p y ) , 
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where 



Fo(p x ,Py) 
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is the electronic Fermi function, and po(t) = — (eEo/fl)e at sm fit is the solution of the single particle equation of 
motion. The electric current j(t) = —eg s g v S v/ p (i) then assumes the form 
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j y = 0, where g s — g v — 2 are the spin and valley degeneracies in graphene, and S is the sample area. After some 
lengthy but simple transformation, Eq. ^ can be rewritten as 



jx{t) = en s V 
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where 



g s gvP% 
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is the density of electrons, pf = £f/V is the Fermi momentum, and 

Q F (t) = = ^^sin(fit) = Q FO sin(Ot) 
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is the field parameter, proportional to the ac electric field Eq. 

Figure [T^, shows the current © as a function of time fit. One sees that in the low-field limit the response is linear. 
Expanding the current ©, we get at Qfo *C 1 



j x (t) S3 en s V<3 
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so that the linear response conductivity (in the collisionlcss approximation) is 
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The expression (1101) co incides with the intra-band Drude conductivity, which can be obtained from the linear-response 
theory [H, [l7|, HI ij. As the inter-band conductivity is of order of e 2 /h 0, E3, S EH, [H [1H , our quasi- 

classical approach is valid at hfl < ef- At the electron density ~ 10 11 — 10 12 cm~ 2 this restricts the frequency by the 
value of 10 - 30 THz. 

In the strong-field limit Qfo 1 Eq. © results in the formula |T]). From the condition Qfo ^ 1, rewritten as 
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FIG. 1: (Color online) The time dependence of the ac electric current, measured in units en s V, at harmonic excitation of the 
system at the frequency £1. (a) The temperature is zero, T/tF = 0; the curves are labeled by the values of the electric field 
parameter Qfo = eEoV/^leF- (b) The temperature is finite, the Fermi energy is zero, ef = 0; the curves are labeled by the 
values of the parameter Qto = eEoV/QT. 



one sees that the required ac electric field grows linearly with the electromagnetic wave frequency and is proportional 
to the square root of the electron density. At / ~ 50 GHz and n s ~ 10 11 cm" 2 , the inequality (JTTJ) is fulfilled at 
Eq ^ 100 V/cm. This value can be reduced in systems with lower electron/hole density. Therefore, we consider now 
an opposite limiting case with ep = 0, but finite temperature T. 

At finite T and the vanishing ep = both electrons and holes contribute to the charge carrier density 



n e + n h = ln o A9 (12) 
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and to the current. Starting again from Eq. ([5]) but accounting for the hole contribution and putting ep = 0, we get 
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where 



QT(t) = -^^ = ^^smnt = Q TO smnt. (14) 

Figure [TJa shows the current (fTB")) as a function of time ftt. In the low-field limit Qto < 1 we get from (fTB"]) the 
current 

j x (t) w en s VQ T (t)^^, (15) 



and the correct expression for the linear-response intra-band dynamic conductivity 17 1, 

One sees that the quasi-classical approach is now valid at Ml < T. This restricts the frequency by the value of ~ 200 
GHz at T ~ 10 K and ~ 6 THz at room temperature. In the strong field regime Qto 1 Eq. (|13p is reduced, again, 
to fTJ). Figure shows the Fourier components of the ac electric current, for m = 1, 3 and 5, as a function of the field 
parameter eEoV/ClT at ep = 0. The strong-field condition now assumes the form 



(17) 
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FIG. 2: (Color online) The Fourier components of the current (|13[) . in arbitrary units, as a function of Qto = eEoV/QT at 
£f/T = 0. 

At T ~ 10 K and / ~ 100 GHz this gives a moderate value of the required electric field E$ ~ 5 V/cm. The efficiency 
of the predicted frequency multiplication effect can be increased further by using the resonance response of the system 
at the plasmon, the cyclotron, or the magnetoplasmon frequency. 

To summarize, we have investigated the non-linear electrodynamic response of 2D electrons and holes in graphene. 
We have shown that irradiation of graphene by the electromagnetic wave with the frequency fi should lead to the 
higher harmonics generation at frequencies 30, 50, et cetera. The efficiency of the frequency up-conversion is rather 
high: the amplitudes of the higher harmonics of the ac electric current fall down slowly (as 1 /m) with the harmonics 
index m. The presented quasi-classical theory is valid at frfl < max{ep,T}. This estimate shows that the effect 
works at frequencies up to 5 — 10 THz, which opens up exciting opportunities for building new graphene devices for 
terahertz and sub-terahertz electronics. 

I wish to thank Klaus Ziegler for stimulating discussion. The work was partly supported by the Swedish Research 
Council and INTAS. 




[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. 

Firsov, Nature 438, 197 (2005). 
[2] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 438, 201 (2005). 
[3] P. R. Wallace, Phys. Rev. 71, 622 (1947). 
[4] G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

[5] K. S. Novoselov, Z. Jiang, Y. Zhang, S. V. Morozov, H. L. Stormer, U. Zeitler, J. C. Maan, G. S. Boebinger, P. Kim, and 
A. K. Geim, Science online (February 15, 2007), cond-mat/0702408. 

[6] V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 146801 (2005). 

[7] K. Ziegler, Phys. Rev. Lett. 97, 266802 (2006). 

[8] M. I. Katsnelson, Europ. Phys. J. B 51, 157 (2006). 

[9] V. V. Cheianov and V. I. Fal'ko, Phys. Rev. Lett. 97, 226801 (2006). 
[10] J. Nilsson, A. H. Castro Neto, F. Guinea, and N. M. R. Peres, Phys. Rev. Lett. 97, 266801 (2006). 
[11] T. Y. Tudorovskiy and A. V. Chaplik, JETP Letters 84, 619 (2006). 
[12] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, 076602 (2007). 
[13] J.-N. Fuchs and P. Lederer, Phys. Rev. Lett. 98, 016803 (2007). 

[14] A. Bostwick, T. Ohta, T. Seyller, K. Horn, and E. Rotenberg, Nature Physics 3, 36 (2007). 

[15] R. S. Deacon, K.-C. Chuang, R. J. Nicholas, K. S. Novoselov, and A. K. Geim (2007), arXiv:0704.0410. 

[16] V. P. Gusynin and S. G. Sharapov, Phys. Rev. B 73, 245411 (2006). 

[17] L. A. Falkovsky and A. A. Varlamov (2006), cond-mat/0606800. 

[18] M. L. Sadowski, G. Martinez, M. Potemski, C. Berger, and W. A. de Heer, Phys. Rev. Lett. 97, 266405 (2006). 
[19] M. L. Sadowski, G. Martinez, M. Potemski, C. Berger, and W. A. de Heer (2007), arXiv:0704.0585. 



■5 



[20] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys. Rev. Lett. 96, 256802 (2006). 

[21] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, J. Phys. Condens. Matter 19, 026222 (2007). 

[22] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte (2007), arXiv:cond-mat/0701053. 

[23] E. H. Hwang and S. Das Sarma (2006), cond-mat/0610561. 

[24] O. Vafek, Phys. Rev. Lett. 97, 266406 (2006). 

[25] V. Apalkov, X.-F. Wang, and T. Chakraborty (2006), cond-mat/0611465. 
[26] D. S. L. Abergel and V. I. Fal'ko (2006), cond-mat/0610673. 
[27] V. Ryzhii, Jpn. J. Appl. Phys. 45, L923 (2006). 

[28] B. Trauzettel, Y. M. Blanter, and A. F. Morpurgo, Phys. Rev. B 75, 035305 (2007). 

[29] V. Ryzhii, A. Satou, and T. Otsuji, J. Appl. Phys. 101, 024509 (2007). 

[30] F. Rana and F. R. Ahmad (2007), arXiv:0704.0607vl [cond-mat.mes-hall]. 

[31] S. A. Mikhailov and K. Ziegler (2007), arXiv:cond-mat/0703406. 

[32] M. I. Katsnelson, Materials Today 10, 20 (2007). 

[33] A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 (2007). 



